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Imposing that the excursion distance of inflaton in field space during inflation be less than the 
Planck scale, we derive an upper bound on the tensor-to-scalar ratio at the CMB scales, i.e. rt,max, in 
the general canonical single-field slow-roll inflation model, in particular the model with non-negligible 
running of the spectral index as and/or the running of running /3s. We find that r,^max — 7 x 10“'* 
for Us = 0.9645 without running and running of running, and r^^rnax is significantly relaxed to the 
order of 0(10”^ ~ 10“^) in the inflation model with Qs and/or /3s ~ -1-0(10“^) which are marginally 
preferred by the Planck 2015 data. 
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Inflation m is the leading paradigm for the early 
universe. Not only does it solve several major puzzles of 
the hot big bang model, such as the flatness puzzle, the 
horizon puzzle and so on, it offers a natural mechanism to 
explain the origin of cosmic structure. The gravitational 
waves can be produced in the early universe [5], whose 
direct detection would be taken as the evidence for the 
inflationary universe and fix the energy scale of inflation. 

In March of 2014, BICEP2 detected an excess of B- 
mode power over the lensed-ACDM expectation on large 
scales at more than 5a CL [6] . Unfortunately, this excess 
does not provide evidence for the detection of the primor¬ 
dial gravitational wave once the contribution to the CMB 
B-mode power from polarized thermal dust was carefully 
estimated in [71H] and Planck HFI 353 GHz CMB polar¬ 
ization data (3 were combined with Planck2013 TT and 
WMAP Polarization data in [10]. Recently the cross¬ 
spectra between BICEP2/Keck Array maps and all the 
polarized bands of Planck (BKP) confirmed that the ex¬ 
cess of B-modes detected by BICEP2 mainly originates 
from the polarized dust m, not the primordial gravi¬ 
tational waves. The combination of the BKP likelihood 
and Planck 2015 likelihood yields the constraint on the 
tensor-to-scalar ratio r: 

’’0.002 < 0.08 (1) 

at 95% CL in [T2|, where the subscript 0.002 means that 
the value of tensor-to-scalar ratio is evaluated at the pivot 
scale 0.002 Mpc“^. 

As an alternative to the inflationary universe, ekpy- 
rotic scenario proposed by Khoury et al. in |13j is also 
supposed to address the flatness, horizon and monopole 
puzzles and generate a nearly scale-invariant spectrum 
of density perturbations. But the ekpyrotic scenario pre¬ 
dicts a strongly blue-tilted tensor power spectrum with 
the spectral index n* = dinPt/dink ~ 2 which pro¬ 
vides an observational signature for distinguishing the 
ekpyrotic scenario from inflation model which predicts 
n* = —r/8 ~ 0 [H], where Pt is the amplitude of the ten¬ 
sor power spectrum. Even though the current data are 


not good enough for tightly constraining the tilt of ten¬ 
sor power spectrum, a strongly blue-tilted tensor power 
spectrum with rit — 2 predicted by the ekpyrotic scenario 
is disfavored at high confidence level m once the con¬ 
straint on the intensity of a stochastic gravitational-wave 
background from the Laser Interferometer Gravitational- 
Waves Observatory (LIGO) [16] was combined. This pro¬ 
vides a strong support to the inflationary universe. 

Since the primordial gravitational waves encode crit¬ 
ical information about inflationary universe, we wonder 
what the predictions about the tensor power spectrum 
in the inflation model are. In Lyth proposed that 
the amplitude of tensor power spectrum is bounded from 
above by imposing a sub-Planckian excursion of inflaton 
in field space during inflation. This is the so-called “Lyth 
bound”. In this letter we will extend the discussion in 
m to more general cases. 

Actually there is a long-term debate about whether 
the excursion distance of inflaton in field space can go 
beyond the Planck scale. 

• The chaotic inflation model with a monomial poten¬ 
tial [18] is a typical large-held inhation model in which 
(j)/Mp ^ 0(10) during inhation. However, according to 
the weak gravity conjecture for the scalar held proposed 
in |19j . the vacuum expectation value of inhaton held 
should be smaller than the Planck scale, and hence we 
concluded that such a model is not self-consistent and 
should be ruled out [19]. Now it has been disfavoured at 
around 3cr CL 

• Another typical large-held inhation model is the natu¬ 
ral inhation model [50], where the inhaton is a pseudo- 
Goldstone boson with a cosine potential. However it re¬ 
quires a super-Planckian decay constant which is not con¬ 
sistent with string theory [5T]. Now this model is also in 
trouble [15] . 

• Brane inhation [22j provides a generic prototype of in¬ 
hation in string theory. In [53] the authors explicitly 
derived the microscopic bound on the maximal held vari¬ 
ation of inhaton during brane inhation in the warped ge¬ 
ometry, and found that it is impossible to achieve a super- 
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Planckian excursion distance of inflaton in field space. 

In this letter we suppose that the excursion distance 
of inflaton in field space during inflation be less than 
the reduced Planck scale Mp, and then derive an upper 
bound on the tensor-to-scalar ratio in the general canon¬ 
ical single-held slow-roll inhation model. 

The dynamics of the canonical single-held slow-roll in¬ 
hation model is govern by 

^ + 3H<p + V'{cl)) = 0 , 

and 


( 2 ) 

(3) 


up to higher order slow-roll parameters, where q ~ 1.063, 
and 




e3 = M, 


VH^) 


£4 = M!^ 


P3(<^) 


-■ (13) 


And the tensor-to-scalar ratio is given by 


r=^- 16ei [1 - 29(2ei - € 2 )] ■ (14) 

See PH and the appendix of PH in detail. Re¬ 
cently Planck collaboration released their scientihc re¬ 
sults about inhation in [12] . The combination of Planck 
TT, TE, EE and lowP datasets implies 


where V{(j)) is the potential of inhaton held and the 
dot and prime denote the derivatives with respect to the 
cosmic time t and the inhaton held (/> respectively. The 
inhaton held slowly rolls down its potential ii (j)^ V 
and \(j)\ ^ or equivalently ei ^ 1 and |e 2 | <C 3, 

where ei and £2 are the slow-roll parameters dehned by 


_ (y'V 

2 \v) ’ 

(4) 

V" 

£2 ^ M^ — . 

(5) 


Us = 0.9645 ± 0.0049, (15) 

at 68% CL. Adding the running of spectral index, they 
found 

=-0.0057 ± 0.0071, (16) 

at 68% CL. However, once the running of running is al¬ 
lowed to hoat, they got 

Us = 0.9586 ±0.0056, (17) 

Os = 0.009 ±0.010, (18) 

f3s = 0.025 ±0.013, (19) 


Usually the power spectrum of the scalar perturbation 
generated by inhation model can be parameterized as 
follows 


fs(fc) — [ -— 




( 6 ) 


where ris,*, as.» and /3s,* are the amplitude of scalar 
power spectrum, the spectral index, the running of spec¬ 
tral index and the running of running at the pivot scale 
fc*. In general, the spectral index of scalar power spec¬ 
trum, the running of spectral index and the running of 
running are dehned by 


ris 

as 


dlnPsjk) 

d\nk 

dus 

dink' 

das 

dink 


(7) 

( 8 ) 
(9) 


at 68% CL. All of the above constraints are evaluated at 
the pivot scale fc* = 0.05 Mpc“^. Even though the scale- 
independent spectral index can ht the data, allowing for 
the running of running provides a better ht to the Planck 
data, such that ~ —4.9. 

In the slow-roll paradigm, the excursion distance of 
inhaton in held space compared to the reduced Planck 
scale is related to the tensor-to-scalar ratio by 

where t at is the cosmic time corresponding to the number 
of e-folds N before the end of inhation, which is dehned 
by A = Hdt. In general, the tensor-to-scalar ratio r 
is a function of time (or equivalently, the e-folding num¬ 
ber before the end of inhation N), and then the bound 
on r for |A(/)|/Mp < I is model-dependent. In [5^, the 
slow-roll parameter ei is supposed to be parametrized by 


In the canonical single-held slow-roll inhation model, 
these three parameters are related to the slow-roll pa¬ 
rameters by 

ris = 1 — 6ei ± 2e2 ± 2qe3 ± • • • , (10) 

as = —246^ ± IheiCi — 2e3 — 2ge4 ± • • • , (H) 

fds = -192e? ± 192e?e2 - 32eie^ - 24eie3 

±2e2e3 ± 2e4 ± • • • , (12) 


ei(A) - eo ± ^2-2£’ (27) 

where cq, c and e are constants, and e G [0,1]. Even 
though this parametrization cannot cover all of the 
single-held inhation model, it is quite generic. For ex¬ 
ample, for eo = 0, = p/2 and £ = 1/2, it corresponds 

to the chaotic inhation model with V{(p) cx eg = 0, 
e = 1/d and c <C 1, it corresponds to the brane inhation 












3 


model with V[(j)) = Vb [l — ^] ■ From Eq. (211, 

requiring \A(j)\/Mp < 1 yields ei(iV) < 1/(2A^^) and then 


r < 


8 

iV2‘ 


( 22 ) 


The CMB scales roughly corresponds to ~ 60 
and then 


r* < 2 X 10"^ 


(23) 


Furthermore, based on this parameterization, we can also 
get the bounds on the spectral index, the running of spec¬ 
tral index and the running of running, namely 




2 


< a. < 0, 


4 

W 


< 13, < 0.(24) 


See [55] in detail. All of these predictions are nicely con¬ 
sistent with the constraints from Planck 2015, for exam¬ 
ple those given in Eqs. (15) and (16). 

Even though a power-law scalar power spectrum can fit 
the Planck data well, there is still a big room for allowing 
a running spectral index, e.g. Eqs. ([^, ([^and([^. In 
order to include these results, we need to go beyond the 
parametrization in Eq. (21). For simplicity, we introduce 
a new variable defined by 


a: = In 


a 

a(t*) ’ 


(25) 


where denotes the time of horizon exit of the pertur¬ 
bation mode fc* during inflation. Since r = Pt/P,, 


dlnr , 

—^ =nt- {n, - 1). (26) 

Considering the consistency relation nt = —r/8, we have 


dlnr 

dx 


K “ 1) + ^ 


(27) 


For ns,» = 0.9645 and r» < 0.08, r increases with the 
expansion of the universe. However, if there are non- 
negligibly positive running of the spectral index and/or 
running of running, n, increases and the right hand side 
of the above equation may flip the sign after a few number 
of e-folds, and then r drops down. Therefore one may 
expect that the Lyth bound is expected to be significantly 
relaxed in the case with non-negligibly positive running 
and/or running or running. Since the perturbation mode 
k exits horizon when k = aH and H is roughly a constant 
during inflation, din A: ~ din a = dx. Integrating over 
Eqs. ([^ and , we obtain 


O^g T 

Tig - Xig :lf (y.g^:^X -\~ ■ 

Therefore 


dlnr 

dx 



(28) 

(29) 


(30) 


Here r is replaced by r* on the right hand side of the 
above equation because |dlnr/da;| <C 1 and r = r* -I- 
r* (dlnr/dx) |r=T-,a^-|-- • • ~ r*. Integrating over the above 
equation, we obtain 


r ~ r*e 


— R{ns,* ,r^ ,0's,* ,0s, 


(31) 


where 




(32) 


(^S,* ~ "g" 


,1 2 I Idg * 3 

X + -Ug^^X P —^x . 

2 b 


Now Eq. (20) reads 


Mr, Vs In 


2 , /^s,* j x') 


(33) 


Requiring that the right hand side of the above equa¬ 
tion be less than one yields the Lyth bound for the gen¬ 
eral canonical single-field slow-roll inflation model. Or 
equivalently, the upper bound on the tensor-to-scalar ra¬ 
tio rIf,max satisfies 


dxexp [-\R{ng^^,r^^max,ag^^,l3gy,x)\^ 


2 ■ 


(34) 


This is the key result of this letter. According to Eq. 
one may worry about that the spectral index significantly 
deviates from unity when a; 1 if the running of spectral 
index and/or running of running are large, and then the 
slow-roll approximation adopted in the previous calcula¬ 
tions is not reliable any more. However, if the running 
of spectral index and/or running of running are positive 
and large, the tensor-to-scalar ratio decreases rapidly and 
becomes quite small after a few e-folding numbers stating 


from the time of t*. In this case the integration in (34) 


mainly comes from integrating over these few e-folding 
numbers during which the slow-roll approximation is still 
applicable. So our formula (34) is still reliable for the 


case with non-negligibly positive running of spectral in¬ 
dex and running of running. 

As a double check, we also adopt the Monte-Carlo sim¬ 
ulation to figure out the Lyth bound. Selecting the initial 
slow-roll parameters at the pivot scale which satisfy the 
constraints from observations and numerically solving the 
flow equations of the slow-roll parameters, i.e. 

~ 2ei(2ei—€ 2 ), (35) 

dx 

^ ~ 2(£-l)eie(; - (£-2)e2e(; - e^+i, (36) 

dx 


where £ > 2, we can generate a vast collection of infla¬ 
tion models in which \A(j)\/Mp < 1 and ei(a:) < 1 for 
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a; S [0,-/V*]. As long as the number of models in such 
a collection is large enough, we can figure out the upper 
bound on the tensor-to-scalar ratio. We will compare the 


Monte-Carlo simulation to the analytic formula (341 in 
the following three cases respectively. 

First of all, let’s consider the case with negligible Ofs,* 
and Ps,*, i.e. ~ Ps,* — 0. If ^ 1, Eq. (34) reads 




n -2 


(37) 


In the limit of rig —>■ 1, — 8/IV^ which is same as 

that in Eq. (22). For W = 60 and = 0.9645, we 
have r* < r^ ^ax — 7 x 10“^. ^ The upper bound on 
the tensor-to-scalar ratio for the case with rig^* < 1 is 
smaller than that for the Harrison-Zel’dovich spectrum 
(n.g_* = 1). The upper bound on the tensor-to-scalar 
ratio is illustrated by the red curve in Fig. We can 



points in Fig. From Fig. we see that all of the 
predictions of these models stay below the bound given 
in Eq. (37). It implies that the formula in Eq. (37) is 
reliable. 

Secondly, we consider the model with non-negligible 
running of spectral index, but ps,* — 0. Taking W = 60, 
we numerically solve Eq. (34) and the results show up 
in Fig. where the red-shaded region comes from the 
uncertainty of rig^* in Eq. (15). The upper bound on 



FIG. 2: The upper bound on the tensor-to-scalar ratio r^^rnax 
in the model with non-negligible running of spectral index. 
The red-shaded region comes from the uncertainty of rig,* in 
Eq. (15l. The grey and light grey bands respectively corre¬ 
spond to the constraints on Og,* at la and 2a CL from Planck 
in Eq. (161. The orange points denote the predictions of a vast 
collection of models generated by the Monte-Carlo simulation, 
in which \ A(f>\/Mp < 1 and ei{x) < 1 for x € [0, A*]. 


FIG. 1: The upper bound on the tensor-to-scalar ratio Vt^max 
is illustrated by the red curve. The grey and light grey bands 
respectively correspond to the constraints on rig,* at la and 
2a CL from Planck in Eq. (151. The orange points denote 


the predictions of a vast collection of models generated by the 
Monte-Carlo simulation, in which \ A(f)\/Mp < 1 and 61 ( 2 ;) < 1 
for X € [0, A*]. 


also check Eq. (37) by using the Monte-Carlo simulation. 


Since Qfg ~ ~ 0 in this case, the slow-roll parameters 

can be truncated as ~ 0 for £ > 3. Following the 
method described in the former paragraph and imposing 
I A^l/Mp < 1 and £ 1 ( 2 ::) < 1 for x G [0, A*], we generate a 
vast collection of models which are denoted by the orange 


^ For Us,* = 0.96, rt^max ~ 6 X 10“^. Similar to 1261 . the authors 
considered a specific assumption of ei(A) scaling as a power of 
1/A, and found a more stringent bound on the tensor-to-scalar 
ratio, i.e. r < 2 X 10“® for ng = 0.96, in m- However, in 
this paper, we do not assume any specific parameterization, and 
hence our results are model-independent. 


the tensor-to-scalar ratio is relaxed to be r:^^rnax — 2 x 
10“^, 1.5 X 10“^ once the 1 and 2 a uncertainties of 
the running of spectral index given in Eq. (16) are taken 
into account. Similar to the previous case, the slow-roll 
parameters are truncated as ~ 0 for £ > 4 and a 
vast collection of models denoted by the orange points 
in Fig. are generated by the Monte-Carlo simulation. 
Since almost all of the orange points do not exceed the 
red-shaded region, the analytic formula (34) is reliable for 
the model with non-negligible running of spectral index. 

Third, we consider the case with both non-negligible 
as * and jSg * ■ Our numerical results are illustrated in 
Fig. [3] where the red-shaded region comes from the un¬ 
certainty of rig,* and ag,* in Eqs. (17) and (18). Taking 
the Icr uncertainty of (3s,* in Eq. (19) into account, we 
find r^^max = 0.25 which says that there is still a big 
room for detecting the primordial gravitational waves 
generated during inflation if the running of running is 
in the order of 4-0(10“^) which is marginally preferred 
by Planck 2015 [T2]. Again, the Monte-Carlo simulation 
(e£ ~ 0, for £ > 5) confirms the reliability of our analytic 
formula (34) in this case. 
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FIG. 3: The upper bound on the tensor-to-scalar ratio rt,max 
in the model with non-negligible running of spectral index 
and running of running. The red-shaded region comes from 
the uncertainty of n^,* and Ofs,* in Eqs. ( |17[ ) and ( |18[ ). The 
grey and light grey bands respectively correspond to the con¬ 
straints on Ps at la and 2a CL from Planck in Eq. (191. The 
orange points denote the predictions of a vast collection of 
models generated by the Monte-Carlo simulation, in which 
lApl/Mp < 1 and ei(a;) < 1 for a; £ [0,Af,]. 


To summarize, we derive the upper bound on the 
tensor-to-scalar ratio for the general canonical single-held 
slow-roll inhation model with sub-Planckian excursion 
distance of inhaton in held space. Comparing with the 
Monte-Carlo simulation, we see that the analytic formula 
(341 works quite well for the case with non-negligible run¬ 
ning of spectral index and/or running of running. With¬ 
out adding the running of running, a power-law scalar 
power spectrum give a good ht to the Planck 2015 data. 
For W = 60 and Ug,* = 0.9645, the tensor-to-scalar ra¬ 
tio is bounded from above by 7 x 10“"^ which might be 
marginally detected in the future. Adding the running 
of running provides a better ht to the Planck data, the 
upper bound is relaxed to be order of 0(0.1), and there is 
still a big room for detecting the primordial gravitational 
waves from both the ground-based experiments and satel¬ 
lites in the near future [H]. See more related discussion 
about the bound on the tensor-to-scalar ratio for the in¬ 
hation model with sub-Planckian excursion in [30j . 
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